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1 Introduction 

High energy nucleus-nucleus collisions provide the means of creating nuclear 
matter in conditions of extreme temperature and density pQ[2j[3]. The kinetic 
energy of the incident projectiles would be dissipated in the large volume of 
nuclear matter involved in the reaction. At large energy or baryon density, 
a phase transition is expected from a state of nucleons containing confined 
quarks and gluons to a state of "deconfined" (from their individual nucleons) 
quarks and gluons, in chemical and thermal equilibrium, covering a volume 
that is many units of the confining length scale. This state of nuclear matter 
was originally given the name Quark Gluon Plasma(QGP) [4J, a plasma being 
an ionized gas. However the results at RHIC [2] indicated that instead of 
behaving like a gas of free quarks and gluons, the matter created in heavy 
ion collisions at nucleon-nucleon cm. energy ^/snn = 200 GeV appears 
to be more like a liquid. This matter interacts much more strongly than 
originally expected, as elaborated in peer reviewed articles by the 4 RHIC 
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experiments [HI El El IB] , which inspired the theorists [9] to give it the new 
name "sQGP" (strongly interacting QGP). 

In the terminology of high energy physics, the QGP or sQGP is called a 
"soft" process, related to the QCD confinement scale 

AQ^-O^GeV^-lfni . (1) 

With increasing temperature, T, in analogy to increasing Q 2 , the strong 
coupling constant a s (T) becomes smaller, reducing the binding energy, and 
the string tension, cr(T), becomes smaller, increasing the confining radius, 
effectively screening the potential |10j: 

3 r 3 r p 

where p = p(T) = 1/ro is the Debye screening mass [10]. For r < 1/p a 
quark feels the full color charge, but for r > the quark is free of the 
potential and the string tension, effectively deconfined. 

There has been considerable work over the past three decades in mak- 
ing quantitative predictions for the QGP [2]. The predicted transition tem- 
perature from a state of hadrons to the QGP varies, from T c ~ 150 MeV 
at zero baryon density, to zero temperature at a critical baryon density 
roughly 1 GeV/fm 3 , ~ 6.5 times the normal density of cold nuclear matter 
(po — 0.14 nucleons/fm 3 , p# ~ 930 MeV), where p# is the Baryon chemical 
potential. A typical expected phase diagram of nuclear matter [11] is shown 
in Fig. [TJ Not distinguished on Fig. [1] in the hadronic phase are the liquid 
self-bound ground state of nuclear matter and the gas of free nucleons [12] . 

A nice feature of the search for the QGP is that it requires the inte- 
grated use of many disciplines in Physics: High Energy Particle Physics, Nu- 
clear Physics, Relativistic Mechanics, Quantum Statistical Mechanics, and, 
recently, AdS/CFT string theory [131 E]- From the point of view of an ex- 
perimentalist there are two major questions in this field. The first is how 
to relate the thermodynamical properties (temperature, energy density, en- 
tropy, viscosity ...) of the QGP or hot nuclear matter to properties that can 
be measured in the laboratory. The second question is how the QGP can be 
detected. 

One of the major challenges in this field is to find signatures that are 
unique to the QGP so that this new state of matter can be distinguished 
from the "ordinary physics" of relativistic nuclear collisions. Another more 
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Figure 1: (left) A proposed phase diagram for nuclear matter [TT]: Temper- 
ature, T, vs Baryon Chemical Potential, /i. 

general challenge is to find effects which are specific to A+A collisions, such 
as collective or coherent phenomena, in distinction to cases for which A+A 
collisions can be considered as merely an incoherent superposition of nucleon- 
nucleon collisions [151 EE E] • 

2 Issues in Relativistic Heavy Ion Physics 

2.1 J/ty suppression — the original "gold-plated" QGP 
signature 

Since 1986, the 'gold-plated' signature of deconfinement was thought to be 
J/^ suppression. Matsui and Satz [18] proposed that J/\l/ production in 
A+A collisions will be suppressed by Debye screening of the quark color 
charge in the QGP. The J/^J is produced when two gluons interact to pro- 
duce a c, c pair which then resonates to form the J/\I/. In the plasma the c, c 
interaction is screened so that the c, c go their separate ways and eventually 
pick up other quarks at the periphery to become open charm. 'Anoma- 
lous suppression" of J/^f was found in Pb+Pb collisions at the CERN SpS 
= 17.2 GeV [19J (e.g. see Fig. [18] below) . This is the CERN fixed tar- 
get heavy ion program's main claim to fame: but the situation is complicated 
because J/\l/ are suppressed in p+A collisions [2D] . 

The search for J/\l/ suppression and thermal photon/dilepton radiation 
from the QGP drove the design of the RHIC experiments. My summary of the 
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different views of dilepton resonances in the High Energy [2Tj and Relativistic 
Heavy Ion [TS] Physics communities since the mid 1980 's is shown in Fig. [2J 
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Figure 2: "The road to success": In High Energy Physics (left) a UA1 
measurement of pairs of muons each with > 3 GeV/c shows two Nobel 
prize winning dimuon peaks and one which won the Wolf prize. Success for 
measuring these peaks in RHI physics is shown schematically on the right. 



2.2 Detector issues in A+A compared to p-p collisions 

Another main concern of experimental design in RHI collisions is the huge 
multiplicity in A+A central collisions compared to p-p collisions. A schematic 
drawing of a collision of two relativistic Au nuclei is shown in Fig. [3^. In 
the center of mass system of the nucleus-nucleus collision, the two Lorentz- 
contracted nuclei of radius R approach each other with impact parameter b. 
In the region of overlap, the "participating" nucleons interact with each other, 
while in the non-overlap region, the "spectator" nucleons simply continue on 
their original trajectories and can be measured in Zero Degree Calorimeters 
(ZDC), so that the number of participants can be determined. The degree 
of overlap is called the centrality of the collision, with b ~ 0, being the most 
central and b ~ 2R, the most peripheral. The maximum time of overlap 
is r D = 2R/ r yc where 7 is the Lorentz factor and c is the velocity of light. 
The energy of the inelastic collision is predominantly dissipated by multiple 
particle production, where n ch , the number of charged particles produced, is 
directly proportional [8] to the number of participating nucleons (Ap art ) as 
sketched on Fig. [3^l. Thus, n c h or the total transverse energy Et in central 
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Figure 3: a) (left) Schematic of collision of two nuclei with radius R and 
impact parameter b. The curve with the ordinate labeled da/dn c h represents 
the relative probability of charged particle multiplicity n c h which is directly 
proportional to the number of participating nucleons, N paTt . b) (right) Trans- 
verse energy (Et) distribution in Au+Au and p-p collisions at y/s^N = 200 
GeV from PHENIX [22]. 

Au+Au collisions is roughly A times larger than in a p-p collision, as shown 
in the measured transverse energy spectrum in the PHENIX detector for 
Au+Au compared to p-p (Fig. [3b) and in actual events from the STAR and 
PHENIX detectors at RHIC in Fig. H 




Figure 4: a) (left) A p-p collision in the STAR detector viewed along the 
collision axis; b) (center) Au+Au central collision at ^/snn = 200 GeV in 
the STAR detector; c) (right) Au+Au central collision at ^snn = 200 GeV 
in the PHENIX detector. 

As it is a daunting task to reconstruct all the particles produced in such 
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events, the initial detectors at RHIC (23] concentrated on the measurement of 
single-particle or multi-particle inclusive variables to analyze RHI collisions, 
with inspiration from the CERN ISR which emphasized those techniques be- 
fore the era of jet reconstruction. There are two major detectors in operation 
at RHIC, STAR and PHENIX, and there were also two smaller detectors, 
BRAHMS and PHOBOS, which have completed their program. As may be 
surmised from Fig. HJ STAR, which emphasizes hadron physics, is most like 
a conventional general purpose collider detector, a TPC to detect all charged 
particles over the full azimuth (A<f> = 2ti) and ±1 units of pseudo-rapidity 
(77), while PHENIX is a very high granularity high resolution special pur- 
pose detector covering a smaller solid angle at mid-rapidity, together with a 
muon-detector at forward rapidity [24J. PHENIX is designed to measure and 
trigger on rare processes involving leptons, photons and identified hadrons at 
the highest luminosities with the special features: i) a minimum of material 
(0.4% X ) in the aperture to avoid photon conversions; ii) possibility of zero 
magnetic field on axis to prevent de-correlation of e + e~ pairs from photon 
conversions; iii) Electro-Magnetic Calorimeter (EMCal) and Ring Imaging 
Cherenkov Counter (RICH) for e ± identification and level- 1 trigger; iv) a 
finely segmented EMCal (£77, 5<p = 0.01 x 0.01) to avoid overlapping show- 
ers due to the high multiplicity and for separation of single-7 and n° up to 
Pt ~ 25 GeV/c; v) EMCal and precison Time of Flight measurement for 
particle identification. 

In addition to the large multiplicity, there are two other issues in RHI 
physics which are different from p-p physics: i) space-time issues, both in mo- 
mentum space and coordinate space — for instance what is the spatial extent 
of fragmentation? is there a formation time/ distance?; ii) huge azimuthal 
anisotropies of particle production in non-central collisions (colloquially col- 
lective flow) which are interesting in their own right but can be troublesome. 

2.3 Collective Flow 

A distinguishing feature of A+A collisions compared to either p-p or p+A 
collisions is the collective flow observed. This effect is seen over the full range 
of energies studied in heavy ion collisions, from incident kinetic energy of 
lOOA MeV to cm. energy of ^snn = 200 GeV [25]. Collective flow, or simply 
flow, is a collective effect which can not be obtained from a superposition of 
independent N-N collisions. 

Immediately after an A+A collision, the overlap region defined by the 
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Figure 5: (left) Almond shaped overlap zone generated just after an A+A 
collision where the incident nuclei are moving along the ±z axis. The reac- 
tion plane by definition contains the impact parameter vector (along the x 
axis) [26]. (right) Measurements of elliptical-flow (i^) for identified hadrons 
plotted as v-i divided by the number of constituent quarks n q in the hadron 
as a function of (a) PT/n q , (b) KET/n q [2 



nuclear geometry is almond shaped (see Fig[5]) with the shortest axis along the 
impact parameter vector. Due to the reaction plane breaking the <fi symmetry 
of the problem, the semi-inclusive single particle spectrum is modified by an 
expansion in harmonics |28j of the azimuthal angle of the particle with respect 
to the reaction plane, (f) — [29J , where the angle of the reaction plane 
is defined to be along the impact parameter vector, the x axis in Fig. |5j 

Ed 3 N d 3 N d 3 N 

dp 3 prdpTdydcj) 2tt pTdpxdy 

The expansion parameter t>2, called elliptical flow, is predominant at mid- 
rapidity. In general, the fact that flow is observed in final state hadrons 
shows that thermalization is rapid so that hydrodynamics comes into play 
before the spatial anisotropy of the overlap almond dissipates. At this early 
stage hadrons have not formed and it has been proposed that the constituent 
quarks flow |30j , so that the flow should be proportional to the number of con- 
stituent quarks n q , in which case v 2/n q as a function of pr/n q would represent 
the constituent quark flow as a function of constituent quark transverse mo- 
mentum and would be universal. However, in relativistic hydrodynamics, at 
mid-rapidity, the transverse kinetic energy, tut — mo = (^t — l)^o = KEt, 
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rather than px is the relevant variable, and in fact v 2 /n q as a function of 
KE T /n q seems to exhibit nearly perfect scaling [27] (Fig. Eb). 

The fact that the flow persists for px > 1 GeV/ c implies that the viscosity 
is small [31], perhaps as small as a quantum viscosity bound from string 
theory [32], rj/s — l/(4ir) where rj is the shear viscosity and s the entropy 
density per unit volume. This has led to the description of the "sQGP" 
produced at RHIC as "the perfect fluid" [9]. 

2.4 Test of Hydrodynamics at LHC 

A test of hydrodynamics at the LHC concerns the possible increase of the 
anisotropic flow v 2 beyond the 'hydrodynamic limit'. Wit Busza's extrap- 
olation |33j of v 2 to the LHC energy is shown in Fig. [6^,, a factor of 1.6 
increase from RHIC. A previous paper by NA49 [M] which compared v 2 
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Figure 6: a) (left) Busza's extrapolation of v 2 to LHC [33]. b) (right) 'Hydro 
Limit' for v 2 /e vs 'Bjorken multiplicity density', (1/ S)dN ch /dy calculated in 
viscous hydrodynamics for several values of the initial energy density eo and 
rj/s > 1/(4tt) [HE]. 

measurements from AGS and CERN fixed target experiments to RHIC as 
a function of the 'Bjorken multiplicity density', (l/S)dN ch /dy, where S = 
is the overlap area of the collision zone, showed an increase in v 2 /e from 
fixed target energies leading to a "hydro limit" at RHIC, where e is the 
eccentricity of the collision zone. This limit was confirmed in a recent cal- 
culation using viscous relativistic hydrodynamics [35] which showed a clear 
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hydro limit of v^/e = 0.20 (Fig. Eb). This limit is sensitive to the ratio of the 
viscosity /entropy density, the now famous rj/s, but negligibly sensitive to the 
maximum energy density of the collision or to (1/ ' S)dN c \ i /dy. Thus, I assume 
that this calculation would give a hydro limit at the LHC not too different 
from RHIC, Vijz ~ 0.20. Busza's extrapolation of a factor of 1.6 increase 
in t> 2 from RHIC to LHC gives v^js = 0.32 at LHC. In my opinion this is a 
measurement which can be done to high precision on the first day of Pb+Pb 
collisions at the LHC, since it is high rate and needs no p-p comparison data. 
Personally, I wonder what the hydro aficionados would say if both Heinz [35] 
and Busza's predictions were correct? 

3 Measurements in p-p collisions at RHIC 

In addition to being the first heavy ion collider, RHIC is also the first po- 
larized proton collider. Proton-proton collisions are performed with both 
beams either longitudinally or transversely polarized EE]- The bunch- 
by-bunch polarization is arranged so that the spin averaged cross section is 
obtained to high accuracy if polarization information is ignored. The empha- 
sis on precision EM calorimetry allows PHENIX to excel in the measurement 
of reactions producing photons, such as direct-single-photon production, or 
particles which decay to photons, tt° — > 7 + 7, r\ — >■ 7 + 7, etc. 

In order to understand whether an effect observed in A+A collisions 
exhibits a sensitivity to collective effects or to the presence of a medium 
such as the QGP it is important to establish a precise baseline measurement 
in p-p collisions at the same value of nucleon-nucleon cm. energy ^snn- 
PHENIX measurements of the invariant cross section, Ed 3 a/dp 3 , for it and 
direct-single-7 production in p-p collisions at y/s = 200 GeV are shown in 
Fig. [37J and Fig. [7b [38J, respectively. The inset on Fig.[7^ shows that the 
7T° cross section is exponential ~ e _6pT for p? < 2 GeV/c, as originally para- 
materized by Cocconi [39J , which is the region of soft-multiparticle physics. 
For pt > 2 GeV/c the spectrum is a power law which is indicative of the 
hard-scattering of the quark and gluon constituents of the proton. The excel- 
lent agreement of the measurements with theory is rewarding, although not 
surprising, since, after all, the discovery of n° production at large transverse 
momentum at the CERN-ISR proved that the partons of deeply inelastic 
scattering (DIS) interacted strongly with each other [391 140] . 
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Figure 7: a) (left) PHENIX measurement of invariant cross section of 7r° vs. 
Pt at mid-rapidity in p-p collisons at y/s = 200 GeV. (37] . b) (right) PHENIX 
measurement of inclusive direct-single 7 in p-p collisons at y/s = 200 GeV 
(o), together with all previous data compared to the theory. |38j 

3.1 The influence of the CERN-ISR 

The ISR discovery [40J (Fig. |8^) showed that the e~ 6pT dependence at low pr 
breaks to a power law with characteristic y/s dependence for pt > 2 GeV/c, 
which is more evident from the log-log plot of subsequent data [H] (Fig. [8b) 
as a function of xt = 2pT/y/s. This plot exhibits that the cross section for 
hard-processes obeys the scaling law: 

where n e ff(xr, y/s) ~ 4 — 6 gives the form of the force-law between con- 
stituents as later predicted by Quantum Chromodynamics (QCD) with non- 
scaling structure and fragmentation functions and running coupling con- 
stant [12] . The more familiar equation for the constituent reaction a + b — > 



10 




Figure 8: a) (left) CCR (4U] measurement of the invariant cross section of 
7T° vs. px at mid-rapidity in p-p collisons for 5 values of y/s. b) (right) 
Later ISR measurement of invariant cross section of tt° vs. xt = 2pr/ 1 \fs at 
mid-rapidity in p-p collisons for 3 values of y/s [H] 

c + d (e.g. g + q — >■ g + q) at parton-parton center-of-mass (cm.) energy \/§ 
in "leading logarithm" pQCD jl3] is: 



-i£/^)A(* 2 )^V(cosn (5) 



dx\dxid cos 6 1 * dsdyd cos 9* 1x\Xi 

where f a (xi), fb{x2), are parton distribution functions, the differential prob- 
abilities for partons a and 6 to carry momentum fractions X\ and X2 of 
their respective protons (e.g. u{x2)) : and where ^* is the scattering an- 
gle in the parton-parton cm. system. The parton-parton cm. energy 
squared is s = X1X2S, where y/s is the cm. energy of the p-p collision. 
The parton-parton cm. system moves with rapidity y = l/21n(xi/x2) in 
the p-p cm. system and the transverse momentum of a scattered parton is 
Pt — Pr — 2 srn ^*- Only the characteristic subprocess angular distribu- 
tions, £ ab (cos6>*) and the coupling constant, a s (Q 2 ) = 12vr/(25 ln(Q 2 / A2 ))> 
are fundamental predictions of QCD [44, 45J. 
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Subsequent ISR measurements utilizing inclusive single or pairs of hadrons 
established that high p T particles in p-p collisions are produced from states 
with two roughly back-to-back jets which are the result of scattering of con- 
stituents of the nucleons as described by Quantum Chromodynamics (QCD), 
which was developed during the course of those measurements. These tech- 
niques have been used extensively and further developed at RHIC since they 
are the only practical method to study hard-scattering and jet phenomena 
in Au+Au central collisions at RHIC energies. 

The di-jet structure of events triggered by a high px tt°, measured via 
two-particle correlations at the ISR, is shown in Fig [9] [46j H7]. The peaks 



0.3<p T <1.0 0.3<p T <1.0 




Figure 9: CCOR [Ml ST] measurements at yfs = 62.4 GeV. a,b) Distributions 
of azimuthal angle (A</>) of associated charged particles of transverse momen- 
tum px a , with respect to a trigger tt° with p^ t > 7 GeV/c, for 5 intervals of 
PT (a y a) (left-most panel) for A<p = ±7r/2 rad about the trigger particle, and 
b) (middle panel) for A<p = ±7r/2 about ir radians (i.e. directly opposite in 
azimuth) to the trigger. The trigger particle is restricted to \r)\ < 0.4, while 
the associated charged particles are in the range \rj\ < 0.7. c) (right panel) 
xe distributions (see text) corresponding to the data of the center panel. 
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on both the same side (Fig. [Hh) as the trigger 7r° and opposite in azimuth 
(Fig.[Hb) are due to the correlated charged particles from jets. The integrated 
(in A0) yield of the away side-particles as a function of the variable xe = 
—pr a cos(A0)/pr t ~ z a /zt, where z t = PT t /pT t is the fragmentation variable of 
the trigger jet (with pT t ) and z a = PT a /pT a is the fragmentation variable of the 
away jet (with pr a ), was thought in the ISR era to measure the fragmentation 
function of the away jet (Fig. [5b ) but was found at RHIC to be sensitive, 
instead, to the ratio of the transverse momenta of the away-jet to the trigger 
jet, x h =p Ta /pT t @E]. 

The QCD subprocess angular distribution £ afe (cos#*) was also first mea- 
sured with two-particle correlations of 7r° pairs of large invariant mass at the 
CERN-ISR @U EU] (Fig. EU]), in agreement with QCD [H US] at a funda- 
mental level. 
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Figure 10: a) (left 3 panels) CCOR measurement [l9j [50] of polar angular 
distributions of 7r° pairs with net pt < 1 GeV/c at mid-rapidity in p-p 
collisions with y'i = 62.4 GeV for 3 different values of tot invariant mass 
M n7T . b) (rightmost panel) QCD predictions for £ afe (cos#*) for the elastic 
scattering of gg, qg, qq', qq, and qq with a s (Q 2 ) evolution. 



3.2 Other ISR discoveries important at RHIC 

Two other ISR discoveries, direct single-7 production and direct-single e 
production, and one near miss, J/\l/ production, are important components 
of physics at RHIC. 
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Direct single-7 production via the inverse QCD-compton process 
g + q — » 7 + g is an important probe in A+A collisions because the 7 is 
a direct participant in the reaction (at the constituent level), which emerges 
from the medium without interacting and can be measured precisely. The 
cross sections for direct single-7 production at yfs = 62.4 GeV [52] are shown 
in Fig. ITTa . Two-particle azimuthal correlations of charged hadrons with 
neutral mesons (tt°), compared to direct-7 (Fig. HTb). show that direct-7 are 
isolated, with no accompanying same-side particles, while ir° have accompa- 
nying particles since they are fragments of jets from high pt partons. 
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Figure 11: a) (left) Compilation of invariant cross sections of direct-7 produc- 
tion at ISR [52]; (right) azimuthal correlations of neutral mesons and direct-7 
with [52]. 

Direct single-e 1 * 1 at a level of e /tt ~ 10~ 4 for all values of \fs at 
the CERN-ISR were discovered before either the J/^> or open-charm [53] 
(Fig. [l~2l . After the discovery of the J in 1974, it was demonstrated that 
the J/'fy was not the source of the single-e 1 * 1 (Fig. [13]) and two years later, 
when open charm was discovered, it was shown that the direct e ± were due 
to the semi-leptonic decay of charm mesons [54] . Fig. [T3k [55] shows the 
first J/^! at the ISR [55], Fig. [T3b shows the best J/^f measurement at the 
ISR [56j while Fig. [T3b [53] shows that the direct electrons (Fig. [12]) are not 
the result of J/* decay since (p T ) = 1.1 ± 0.05 GeV/c [56]. 
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Figure 13: a) (left) First J/V at ISR [52]; b) (center) Best 
da ee /dm ee dy\ y= o [56]; c) (right) direct-e 1 * 1 data at y/s = 52.7 GeV (Fig. [12]) 
with calculated e ± spectrum for J/ty for several values of (px) [53] . 

4 From ISR p-p to RHIC A+A physics 

Since hard-scattering at high px > 2 GeV/c is point-like, with distance scale 
1/Pt < 0.1 fm, the cross section in p+A (A+A) collisions, compared to p- 
p, should be larger by the relative number of possible point-like encounters, 
a factor of A (A 2 ) for p+A (A+A) minimum bias collisions. When the 
impact parameter or centrality of the collision is defined, the proportionality 
factor becomes (Taa), the average overlap integral of the nuclear thickness 
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functions. 



4.1 Jet quenching from inclusive it production 
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Figure 14: a) (left) Ed 3 a j dp 3 vs. x T for PHENIX mid-rapidity 7r° at 



= 200 GeV in p-p collisions [37] plus PHENIX [57] and CCOR-ISR 
measurements at y/s = 62.4 GeV, where the absolute px scale of the ISR 
measurement has been corrected upwards by 3% to agree with the PHENIX 
data, b) (right) n° p-p data vs. px at yfs = 200 GeV from a) multiplied by 
(Taa) for Au+Au central (0-10%) collisions compared to semi-inclusive 7r° 
invariant yield in Au+Au central (0-10%) collisions at a/s/vat = 200 GeV. 



The discovery, at RHIC, that it are suppressed by roughly a factor of 
5 compared to point-like scaling of hard-scattering in central Au+Au colli- 
sions is arguably the major discovery in Relativistic Heavy Ion Physics. In 
Fig. H4"k). the PHENIX measurement of Ed 3 a/dp 3 for tt° production in p- 
p collisions at yjs = 62.4 GeV [57] is in excellent agreement with the ISR 
data and the PHENIX vr° data follow the same trend as the lower energy 
data, with a pure power law, Ed 3 a j dp 3 oc p^ 8 ' 1±01 for p T > 3 GeV/c at 
y/s = 200 GeV. In Fig. Ob), the 200 GeV p-p data, multiplied by the point- 
like scaling factor (Taa) f° r (0-10%) central Au+Au collisions are compared 
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to the semi-inclusive invariant ir° yield in central (0-10%) Au+Au collisions 
at y/SNN = 200 GeV and, amazingly, the Au+Au data follow the same 
power-law as the p-p data but are suppressed from the point-like scaled p-p 
data by a factor of ~ 5, independent of pr- The suppression is represented 
quantitatively by the "nuclear modification factor" , Raa(pt), the ratio of the 
measured semi-inclusive yield in A+A collisions to the point-like scaled p-p 
cross section at a given p?'- 

d 2 Nl A /dp T dyN A A 



Raa(pi 



{T A A)d 2 aZ Jdp T dy 



(6) 



In Fig. IT5"a . Raa(pt) is shown for n°, r\ mesons and direct-7 for ^/snn = 
200 GeV Au+Au central (0-10%) collisions. The 7r° and r\ mesons, which are 
fragments of jets from outgoing partons are suppressed by the same amount 
while the direct-7 which do not interact in the medium are not suppressed. 
This indicates a strong medium effect on outgoing partons. Fig. [T5b shows 
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Figure 15: a) (left) Nuclear modification factor, Raa for direct-7, 7r° and rj 
in Au+Au (0-10%) central collisions at ^s^n = 200 GeV [58], together with 
GLV theory curve [59] . b) PHENIX Raa for ^° m Cu+Cu central collisions at 
= 200, 62.4 and 22.4 GeV [60], together with Vitev theory curves [HI]. 



that Raa for central (0-10%) Cu+Cu collisions is comparable at ^snn = 62.4 
and 200 GeV, but that there is no suppression, actually a Cronin enhance- 
ment, at ^s NN = 22.4 GeV. This indicates that the medium which sup- 
presses jets is produced somewhere between y/s^N = 22.4 GeV, the SpS 
Fixed Target highest cm. energy, and 62.4 GeV. 

The measurements at RHIC appear to be in excellent agreement with 
the theoretical curves [59, 61J. The suppression can be explained by the 
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energy loss of the outgoing partons in the dense color- charged medium due 
to coherent Landau-Pomeranchuk-Migdal radiation of gluons, predicted in 
QCD [62], which is sensitive properties of the medium. Measurements of two- 
particle correlations (discussed below, Sec. [7J confirm the loss of energy of 
the away-jet relative to the trigger jet in Au+Au central collisions compared 
to p-p collisions. However, lots of details remain to be understood. 



5 Direct photons at RHIC: thermal photons? 

5.1 Internal Conversions — the first measurement any- 
where of direct photons at low 

Internal conversion of a photon from 7r° and rj decay is well-known and is 
called Dalitz decay [63]. Perhaps less well known in the RHI community is 
the fact that for any reaction (e.g. q + g — > 7 + q) in which a real photon can 
be emitted, a virtual photon (e.g. e + e~ pair of mass m ee > 2m e ) can also 
be emitted. This is called internal-conversion and is generally given by the 
Kroll-Wada formula 



1 dN ee 2a 1 m 2 e 3 



iV 7 dm ee 3tt m ee M 2 



x 



\F(m 



\ 



+ , (7) 



GG 



where M is the mass of the decaying meson or the effective mass of the 
emitting system. The dominant terms are on the first line of Eq. [7J the 
characteristic l/m ee dependence; and the cutoff of the spectrum for m ee > M 
(Fig. [T6"a) [65] . Since the main background for direct-single-7 production is 
a photon from 7r° — > 7 + 7, selecting m ee >100 MeV/c 2 effectively reduces the 
background by an order of magnitude by eliminating the background from 
7T° Dalitz decay, tt° — > 7 + e + + e~ , at the expense of a factor ~ 1000 in rate. 
This allows the direct photon measurements to be extended (for the first 
time in both p-p and Au+Au collisions) below the value of px ~ 4 GeV/c, 
possible with real photons, down to pt — 1 GeV/c (Fig. [16b) [65], which is 
a real achievement. The solid lines on the p-p data are QCD calculations 
which work down to pt = 2 GeV/c. The dashed line is a fit of the p-p data 
to the modified power law B(l + p 2 -./b)~ n , used in the related Drell-Yan [66] 
reaction, which flattens as px — > 0. 
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The relatively flat, non-exponential, spectra for the direct-7 and Drell- 
Yan reactions as pt — > is due to the fact that there is no soft-physics 
production process for them, only production via the partonic subprocesses, 
g + q — > 7 + q and q + q — > e + + e~, respectively This is quite distinct from 
the case for hadron production, e.g. tt°, where the spectra are exponential as 
Pt — > in p-p collisions (Fig. [Tfet) due to soft-production processes, as well as 
in Au+Au collisions. Thus, for direct-7 in Au+Au collisions, the exponential 
spectrum of excess photons above the (Taa) extrapolated p-p fit is unique 
and therefore suggestive of a thermal source. 




Figure 16: a) (left) Invariant mass (m ee ) distribution of e + e~ pairs from 
Au+Au minimum bias events for 1.0 < px < 1.5 GeV/c [65]. Dashed lines 
are Eq. [7] for the mesons indicated. Blue solid line is / c (m), the total di- 
electron yield from the sum of contributions or 'cocktail' of meson Dalitz 
decays; Red solid line is fdir{m) the internal conversion m ee spectrum from 
a direct-photon (M >> m ee ). Black solid line is a fit of the data to the sum 
of cocktail plus direct contributions in the range 80 < m ee < 300 MeV/c 2 . 
b) (right) Invariant cross section (p-p) or invariant yield (Au+Au) of direct 
photons as a function of pt [65] . Filled points are from virtual photons, open 
points from real photons. 
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5.2 Low px vs high px direct-7- 
busy plot 



Learn a lot from a 



The unique behavior of direct-7 at low pt in Au+Au relative to p+p com- 
pared to any other particle is more dramatically illustrated by examining 
the Raa of all particles measured by PHENIX in central Au+Au collisions 
200 GeV (Fig. ED [67]. For the entire region p T < 20 GeV/c 



at y/SNN 

so far measured at RHIC, apart from the p + p which are enhanced in the 
region 2 < Pt^,4 GeV/c ('the baryon anomaly'), the production of no other 
particle is enhanced over point-like scaling. The behavior of Raa of the low 
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Figure 17: Nuclear Modification Factor, Raa(pt) for all identified parti- 
cles so far measured by PHENIX in central Au+Au collisions at a/J/wv = 
200 GeV. [67] 

Pt < 2 GeV/c direct-7 is totally and dramatically different from all the 
other particles, exhibiting an order of magnitude exponential enhancement 
as pt 0. This exponential enhancement is certainly suggestive of a new 
production mechanism in central Au+Au collisions different from the con- 
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ventional soft and hard particle production processes in p-p collisions and its 
unique behavior is attributed to thermal photon production by many authors 
(e.g. see citations in reference |65j). 

5.2.1 Direct photons and mesons up to p T = 20 GeV/c 

Other instructive observations can be gleaned from Fig. [T7J The tt° and r] 
continue to track each other to the highest p?- At lower p?, the meson 
tracks the K ± very well, but with a different value of Raa{pt) than the 7r°, 
while at higher p^,the <ft an d w vector mesons appear to track each other. 
Interestingly, the J/^> seems to track the ir° for < pt < 4 GeV/c; and it 
will be important to see whether this trend continues at higher pj>. 

6 J/^f suppression, still golden? 

The dramatic difference in tt° suppression from SpS to RHIC cm. en- 
ergy (Fig. [15b) is not reflected in J/\I/ suppression, which is nearly iden- 
tical at mid-rapidity at RHIC compared to the NA50 measurements at SpS 
(Fig. HSb) [681 [69] . This casts new doubt on the value of J/ty suppression 
as a probe of deconfinement in addition to the previous complication that 
J/ty are already suppressed (compared to point-like scaling) in p+A and 
B+A collisions (Fig. [T8a). One possible explanation is that c and c quarks 
in the QGP recombine to regenerate J/^f, miraculously making the observed 
Raa equal at SpS and RHIC cm. energies (Fig. [T9"a) [69, 70J. The good 
news is that such models predict the vanishing of suppression or even 
an enhancement (Raa > 1) at LHC energies [711 H2J [73], which would be 
spectacular, if observed. 

Even without the LHC startup, there has been progress this past year 
when, after ~ 20 years (!), p+A comparison data for the J/^f from the 
CERN fixed target program at 158A GeV/c finally became available [74] . 
The cold nuclear matter effect of J/vP suppression in p+A collisions is pa- 
rameterized by an effective absorption cross section cr a [ s which had been 
previously measured to be 4.3 ± 1.0 mb at 400 GeV/c proton beam energy 
and "assumed to be independent of beam energy". The actual measurement 
for 158 GeV p+A collisions gives a J J^ = 7.6 ± 0.9 mb which considerably 
reduces the "anomalous suppression" effect shown in Fig. [TSh to such an ex- 
tent that there is now a clear difference between the CERN SpS and RHIC 
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Figure 18: a) (left) Total cross section for J/ty production divided by AB in 
A+B collisions at 158-200^4 GeV [19] . b) (right) J/ty suppression relative to p-p 
collisions (Raa) as a function of centrality (A^ part ) at RHIC [681 EH] an d at the 
CERN/SPS pi]. 
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Figure 19: a) (left) PHENIX measurement of Raa as a function of centrality from 
Fig. [T8b together with prediction from a coalescence model [70] : b) (right) Raa 
of J/\E' at SpS and RHIC cm. energies normalized to the measured Raa(CNM) 
from cold nuclear matter 1751 
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J/\& suppression for the most central A+A collisions relative to the measured 
Cold Nuclear Matter effect (Fig. IT§b) [75] . Maybe there is still some hope 
for J/^ suppression as a QGP signature, but there is an important lesson 
for LHC. Comparison data for p-p and p+A MUST be taken and must be 
at the same ^/snn as the A+A data. 



7 Two-particle correlations 

If the 7T° suppression shown in Fig. [15] is in fact explained by the energy loss 
of the outgoing partons in the dense color-charged medium, this can be con- 
firmed by measurements of two-particle correlations. These measurements 
are sensitive to the ratio of the energy of the away-jet to the trigger jet, which 
can be compared in Au+Au collisions and p-p collisions. In analogy to Fig. [9] 
(above), the two-particle correlations in Au+Au collisions (Fig. 120k ) show 
clear di-jet structure in both peripheral and central collisions. The away- 
side correlation in central Au+Au collisions is much wider than in peripheral 
Au+Au and p-p collisions and is further complicated by the large multipar- 
ticle background which is a modulated in azimuth by the t>2 collective flow of 
a comparable width to the jet correlation. After the t>2 correction, a double 
peak structure ~ ±1 radian from 7r is evident, with a dip at ir radians. This 
may indicate a reaction of the medium to a passing parton in analogy to a 
"sonic-boom" [7B] and is under active study both theoretically and experi- 
mentally. The energy loss of the away-parton is indicated by the fact that the 
xe distribution in Au+Au central collisions (Fig. [20b) is steeper than that 
from p-p collisions. As noted above, we found in PHENIX [4_8j [77J that the 
Xe distribution did not measure the fragmentation function of the away-jet 
but is sensitive instead to Xh, the ratio of the transverse momentum of the 
away-parton to that of the trigger parton, specifically [35] : 

dP 



dxi 



= JV(n- 1)^-7; 1 , (8) 

PTt >x h {i + x E /x h y 



where N is a normalization factor, and n (=8.1 at 200 GeV) is power of the 
inclusive invariant p? t distribution. 
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Figure 20: a) (left) Azimuthal correlation C(A0) of h^ 1 with 1 < px a < 2.5 
GeV/c with respect to a trigger with 2.5 < pr t < 4 GeV/c in Au+Au: 
(top) central collisions, where the line with data points indicates C(A0) 
before correction for the azimuthally modulated (V2) background, and the 
other line is the v 2 correction which is subtracted to give the jet correlation 
function J(A0) (data points); (bottom)-same for peripheral collisions, b) 
(right) xe ~ PT a /pTt distribution for the Au+Au-central data compared to 
p-p. 

8 A charming surprise 

We designed PHENIX specifically to be able to detect charm particles via 
direct-single e since this went along naturally with J/^ — > e + + e~ de- 
tection and since the single particle reaction avoided the huge combinatoric 
background in Au+Au collisions. We thought that the main purpose of 
open charm production, which corresponds to a hard-scale (m cg >3 GeV/c 2 ), 
would be a check of our centrality definition and (Taa) calculation since the 
total production of c quarks should follow point-like scaling. In fact, our 
first measurement supported this beautifully [78J. However, our subsequent 
measurements proved to be much more interesting and even more beauti- 
ful. Figure l2"Ta shows our direct-single-e 111 measurement in p-p collisions at 
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y/s = 200 GeV |79j in agreement with a QCD calculation of c and b quarks 
as the source of the direct-single-e^ 1 (also called non-photonic e ± at RHIC). 
The total yield of direct-e^ 1 for > 0.3 GeV/c was taken as the yield of 




Figure 21: a) (left) Invariant cross section of direct in p-p collisions [79] 
compared to theoretical predictions from c and b quark semileptonic decay, 
b) (right) Raa as a function of centrality (N paiil ) for the total yield of e ± 
from charm (p? > 0.3) GeV/c, compared to the suppression of the e yield 
at large pt > 3.0 GeV/c which is comparable to that of tt° with (jp? > 4 
GeV/c) [7J 

c-quarks in p-p and Au+Au collisions. The result, Raa = 1 as a function of 
centrality (Fig. l2"Tb). showed that the total c— (c) production followed point- 
like scaling, as expected. The big surprise came at large px where we found 
that the yield of direct-single-e 111 for pt > 3 GeV/c was suppressed nearly 
the same as the 7r° from light quark and gluon production. This strongly 
disfavors the QCD energy-loss explanation of jet-quenching because, naively, 
heavy quarks should radiate much less than light quarks and gluons in the 
medium; but opens up a whole range of new possibilities including string 
theory [SO]. 

The suppression of direct-single-e 1 * 1 is even more dramatic as a function 
of pr^,5 GeV/c (Fig 122k) which indicates suppression of heavy quarks as 
large as that for ir° in the region where the m>4 GeV 6-quarks dominate. 
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Figure [22b shows that heavy quarks exhibit collective flow (t^), another 
indication of a very strong interaction with the medium. 
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Figure 22: a) (left) Raa (central Au+Au) b) (right) t>2 (minimum bias 
Au+Au) as a function of pt for direct-e 1 * 1 at a/s/vjv = 200 GeV 



9 Zichichi to the rescue? 

In September 2007, I read an article by Nino, "Yukawa's gold mine" in the 
CERN Courier taken from his talk at the 2007 International Nuclear Physics 
meeting in Tokyo, Japan, in which he proposed: "We know that confinement 
produces masses of the order of a giga-electron-volt. Therefore, according 
to our present understanding, the QCD colourless condition cannot explain 
the heavy quark mass. However, since the origin of the quark masses is still 
not known, it cannot be excluded that in a QCD coloured world, the six 
quarks are all nearly massless and that the colourless condition is 'flavour' 
dependent." 

Nino's idea really excited me even though, or perhaps because, it appeared 
to overturn two of the major tenets of the Standard Model since it seemed to 
imply that: QCD isn't flavor blind; the masses of quarks aren't given by the 
Higgs mechanism. Massless b and c quarks in a color-charged medium would 
be the simplest way to explain the apparent equality of gluon, light quark 
and heavy quark suppression indicated by the equality of Raa for n° and 
direct single-e 1 * 1 in regions where both c and b quarks dominate. Furthermore 
RHIC and LHC-Ions are the only place in the Universe to test this idea. 
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It may seem surprising that I would be so quick to take Nino's idea 
so seriously. This confidence dates from my graduate student days when I 
checked the proceedings of the 12th ICHEP in Dubna, Russia in 1964 to see 
how my thesis results were reported and I found several interesting questions 
and comments by an "A. Zichichi" printed in the proceedings. One comment 
about how to find the W boson in p+p collisions deserves a verbatim quote 
because it was exactly how the W was discovered at CERN 19 years later: 
"We would observe the /x's from W-decays. By measuring the angular and 
momentum distribution at large angles of K and 7r's, we can predict the 
corresponding /z-spectrum. We then see if the /x's found at large angles agree 
with or exceed the expected numbers." 

Nino's idea seems much more reasonable to me than the string theory 
explanations of heavy-quark suppression (especially since they can't explain 
light-quark suppression). Nevertheless, just to be safe, I asked some distin- 
guished theorists what they thought, with these results: 

• Stan Brodsky: "Oh, you mean the Higgs field can't penetrate the QGP." 

• Rob Pisarski: " You mean that the propagation of heavy and light 
quarks through the medium is the same." 

• Chris Quigg (Moriond 2008): "The Higgs coupling to vector bosons 7, 
W, Z is specified in the standard model and is a fundamental issue. 
One big question to be answered by the LHC is whether the Higgs gives 
mass to fermions or only to gauge bosons. The Yukawa couplings to 
fermions are put in by hand and are not required." "What sets fermion 
masses, mixings?" 

• Bill Marciano: "No change in the t-quark, W, Higgs mass relationship 
if there is no Yukawa coupling: but there could be other changes." 

Nino proposed to test his idea by shooting a proton beam through a QGP 
formed in a Pb+Pb collision at the LHC and seeing the proton 'dissolved' by 
the QGP. My idea is to use the new PHENIX VTX detector, to be installed 
in 2010, to map out, on an event-by-event basis, the di-hadron correlations 
from identified b—b di-jets, identified c—c di-jets, which do not originate from 
the vertex, and light quark and gluon di-jets, which originate from the vertex 
and can be measured with 7r°-hadron correlations. A steepening of the slope 
of the xe distribution of heavy-quark correlations as in Fig. [20b will confirm 
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in detail (or falsify) whether the different flavors of quarks behave as if they 
have the same energy loss (hence mass) in a color-charged medium. If Nino's 
proposed effect is true, that the masses of fermions are not given by the Higgs 
particle, and we can confirm the effect at RHIC or LHC-Ions, this would be a 
case where we Relativistic Heavy Ion Physicists may have something unique 
to contribute at the most fundamental level to the Standard Model, which 
would constitute a "transformational discovery." Of course the LHC could 
falsify this idea by finding the Higgs decay to b — b at the expected rate in 
p-p collisions. Clearly, there are exciting years ahead of us! 
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